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Features of particle emission and critical point behavior are investigated as functions of the isospin 
of disassembling sources and temperature at a moderate freeze-out density for medium-size Xe 
isotopes in the framework of isospin dependent lattice gas model. Multiplicities of emitted light 
particles, isotopic and isobaric ratios of light particles show the strong dependence on the isospin 
of the dissociation source, but double ratios of light isotope pairs and the critical temperature 
determined by the extreme values of some critical observables are insensitive to the isospin of the 
systems. Values of the power law parameter of cluster mass distribution, mean multiplicity of 
intermediate mass fragments (lAIF), information entropy (_ff) and Campi's second moment (&) 
also show a minor dependence on the isospin of Xe isotopes at the critical point. In addition, the 
slopes of the average multiplicites of the neutrons (Nn), protons (Np), charged particles (Ncp), and 
IMFs (Nimf), slopes of the largest fragment mass number (Amax), and the excitation energy per 
nucleon of the disassembling source {E* /A) to temperature are investigated as well as variances of 
the distributions of 7V„, Np, Ncp, Nimf, Amax and E* /A. It is found that they can be taken as 
additional judgements to the critical phenomena. 

PACS Number(s): 25.70.Pq, 05.70.Jk, 24.10.Pa, 24.60.Ky 



I. INTRODUCTION 

Isospin influence on the formation and decay of hot nu- 
clei is an important subject in heavy ion collision physics 
nowadays. Interests in this direction have been largely 
pushed away with the development of the radioactive 
beam technique. Many new impressive experiments aim- 
ing to explore such isospin effects can be performed by us- 
ing the radioactive beams and/or targets with large neu- 
tron or proton excess. It offers the possibility to study the 
properties of nuclear matter in the range from symmet- 
rical nuclear matter to pure neutron matter. Recently 
some theoretical investigations to the nuclear equation 
of state, chemical and mechanical instabilites as well 
as liquid-gas phase transition for isospin asymmetrical 
nuclear matter were performed already ||l|. In addi- 
tion, the study on isospin dependent nucleon - nucleon 
cross section [^-^ is also an important subject due to 
its significant effect on the dynamical process of heavy 
ion reactions induced by radioactive beams. Experimen- 
tally some new isospin dependent phenomena have been 
also discovered. For examples, the iso&™n dependences 
of preequilibrium nucleon emission] nuclear stop- 
ping ||l^, nuclear collective flow ]ri,|l2|, total reaction 
cross section, radii of neutron-rich nuclei |p^-[l5|| and sub- 
threshold pion production |l6| have been studied by sev- 
eral groups. However, more experimental and theoretical 
studies are still needed to better understand the isospin 



physics. 

On the other hand, phase transition and critical point 
behavior in a finite nuclear system is a debating topics. 
Two kind of important phase transitions of nuclear mat- 
ter, namely the liquid gas phase transition and the quark- 
gulon plasma phase transition which have been predicted 
to occur in HIC at intermediate and ultra-relativistic en- 
ergy domain, respectively, are attracting more and more 
nuclear scientists to pay attentions. Due to the man- 
ifest significance of both phase transitions in clarifing 
the properties of nuclear matter in extreme conditions, 
theoretical physicists are trying to present various possi- 
ble judgements to characterize the phase transition with 
divers models as well as experimentalists are searching 
for such any evidences with the help of advanced and 
complicated set-up. 

In this article, we will investigate some isospin effects 
on cluster emission and critical point behavior stemming 
from nuclear liquid gas phase transition in the framework 
of isospin dependent lattice gas model. The paper is or- 
ganized as following: In Sec. II, a brief description of 
the isospin dependent lattice gas model is given. Results 
and discussions for the disassembly of I22,i29,i37,i46^g 
isotopes in a fixed freeze-out density of ^ 0.39po are 
presented in Sec. III. Firstly, the multiplicities of particle 
emission are used to discuss the influence of the isospin. 
In addition, the slopes of some average quantities to tem- 
perature and the variances of the distributions of some 
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physical observables are explored to be the judgements 
to locate critical point. Secondly, the isospin dependence 
of the ratios between two light particles (LP) is investi- 
gated. Thirdly, the double ratios between two pair light 
isotopes are studied. Fourthly, the critical observables 
with different isospin are presented. Finally the sum- 
mary is given in Sec. IV. 



II. DESCRIPTION OF MODEL 

The lattice gas model of Lee and Yang , where the 
grancanonical partition function of a gas with one type of 
atoms is mapped into the canonical ensemble of an Ising 
model for spin 1/2 particles, has successfully described 
the liquid-gas phase transition for atomic system. The 
same model has already been applied to nuclear physics 
for isospin symmetrical systems in the grancanonical en- 
semble fl^ with an approximate sampling of the canon- 
ical ensemble [p"9|-p3t, and also for isospin asymmetrical 
nuclear matter in the mean field approximation p^ . In 
this article, we will adopt the similar lattice gas model 
developed by J. Pan and S. Das Gupta 21|. Some 
details and features of the model can be found in their 
papers. In comparison with the earlier version of their 
model 1^,^, the different points in this work are that 
the exact Metropolis sampling for placing nucleons on 
the cubic lattice are adopted and while the isospin de- 
pendent interaction potential is taken (which has been 
also incarnated in their recent works [^^). To better 
understand the context of this work we will make a brief 
description for the model and Monte-Carlo Metrolpois 
sampling technique. 

In the lattice gas model, A nucleons with an occupa- 
tion number s which is defined s = 1 (-1) for a proton 
(neutron) or s = for a vacancy, are placed on the L 
sites of lattice. Nucleons in the nearest neighboring sites 
have interaction with an energy e^^s • The hamiltonian 
is written by 



i=l 



2m 



SiSj 



(1) 



1<J 



The interaction constant e^.g . is fixed to reproduce the 
binding energy of the nuclei. In order to incorporate 
the isospin effect in the lattice gas model, the short- 
range interaction constant tsiSj is made the difference be- 
tween the nearest neighboring like-nucleons and unlike- 



nucleons: 



en = 0. MeV, 



Enp = £1-1 = £-11 — - 5.33 MeV, which indicates the 
repulsion between the nearest neighboring like-nucleons 
and attraction between the nearest neighboring unlike- 
nucleons. This kind of isospin dependent interaction in- 
corporates, in a certain extent, Pauli exclusion principal 
and avoids effectively to produce some unreasonable clus- 
ters, like di-proton and di-neutron etc. Three-dimension 
cubic lattice with L sites is used which results that an 



assumed freeze-out density of disassembling system is p/ 
= ^Poj in which po is the normal nuclear density. The 
disassembly of the system is to be calculated at p/, be- 
yond which nucleons are too far apart to interact. N + Z 
nucleons are put into L sites by Monte Carlo sampling 
using the exact Metropolis algorithm [|9|j2^ . 

As pointed out in Ref. |^,^, one has to be careful 
to treat the process of Metropolis sampling to fulfill the 
detailed balance and therefore warrant the correct equi- 
librium distribution in final state. Speaking in detail, in 
this work, first we establish an initial configuration with 
N + Z nucleons. Second, for each event, we will test suf- 
ficient number of "spin" -exchange steps, eg. 20000 steps 
in this work to let finally the system generates states with 
a probability proportional to the Boltzmann probability 
distribution with Metropolis algorithm. In each "spin"- 
exchange step, we make a random trial change on the 
basis of the previous configuration. For example choose a 
nucleon at random and attempt to exchange it with one 
of its neighboring nucleons or the vacancies at random 
regardless of the sign of its "spin" (Kawasaki- like spin- 
exchange dynamics |3^), then compute the change /S.E 
in the energy of the system due to the trial change. If /S.E 
is less than or equal to zero, accept the new configuration 
and repeat the next "spin" -exchange step. If Ai? is posi- 



tive, compute the "transition probability" W 



and compare it with the a random number r in the in- 
terval [0,1]. If r < W , accept the new configuration; oth- 
erwise retain the previous configuration, and then start 
the next "spin" -exchange step. 20000 " spin" -exchange 
steps are performed to assure to get the equilibrium state. 
Third, once the nucleons have been placed stably on the 
cubic lattice after 20000 "spin" -exchange steps for each 
event, their momenta are generated by Monte Carlo sam- 
pling of Maxwell Boltzmann distribution. Thus various 
observables can be calculated in a straightforward fashion 
for each event. One point should be emphasized here is 
that the above Monte-Carlo Metropolis " spin" -exchange 
approach between the nearest neighbors, independent of 
their "spin", is evidenced to be satisfied to the detailed 
balance as noted in Ref. ||2^,|2^ . In other word this sam- 
pling method will guarantee that the generated micro- 
scopic states form an equilibrium canonical ensemble. 

One of the basic measurable quantities is the size dis- 
tribution of clusters. The definition of clusters has been 
extensively discussed in the previous references [p0|-p2|. 
The first method is so-called Ising cluster which com- 
bines all the connected sites as cluster. However, it is 
not proper approach to define clusters in the lattice gas 
model since it does not fulfill the requirement that the 
correlation length should diverge at the critical point . 
In condensed matter physics, Coniglio and Klein ||3l[] pro- 
posed to combine the above site percolation with an ad- 
dition bond percolation algorithm using a temperature 
dependent bonding probability p(r) = ^_ ^-<^!!iajS,sj /2T ^ 
In the lattice gas model a similar way to Coniglio-Klein's 
prescription is extensively adopted to define the clus- 
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ter [po|-p3[: ie. two neighboring nucleons are viewed 
to be in the same fragment if their relative kinetic en- 
ergy is insufficient to overcome the attractive potential: 
— CsiSjSiSj < 0. ft results in a similar tempera- 
ture dependent bonding probability to Coniglio-Klein's 
prescription. 

III. RESULTS AND DISCUSSIONS 

Four isotopes of Xe are chosed as examples to illustrate 
the isospin effect with the help of the isospin dependent 
lattice gas model. Their isospin parameter (/ = ^^^ ) is 
about 0.11, 0.16, 0.21 and 0.26 for ^^^Xe, ^^^Xe, "7Xe 
and -^^^Xe, respectively. The freeze-out density pf has 
been chosen to be close to 0.39 po, as extracted from the 
analysis of Ar Sc |§ and ^^Cl + Au and ™Ge + Ti 
p2[ with the same model. There is also good support 
from experiment that the value of p/ is significantly be- 
low 0.5/90 [Q- 343 sites with 7x7x7 cube results in 
0.36, 0.38, 0.40, and 0.43 po of the freeze-out density pf 
fQi- I22,i29,i37,i46j^g^ respectively. Calculations were per- 
formed from 3 to 7 MeV with a step 0.5 McV and 1000 
events were accumulated at each temperature for each 
isotope. 

A. Cluster Emission and Its Application to Locate 
the Critical Temperature 

1. Multiplicity of Clusters, the Largest Fragment Mass and 
Excitation Energy 

Before presenting the calculated results, we firstly in- 
troduce a definition of excitation energy in this lattice 
gas model. The excitation energy per nucleon can be 
written as |^ 

3 iV^ Ns-"- 

E*/A = ET-E,.^.^{-T + e,,,^)-e^,^ (2) 

in which N^^^^' and N^^ is the number of the bonds of 
unlike-nucleons in the ground state (zero temperature) 
and at T, respectively. Experimentally the excitation 
energy of a nuclear system is generally given with re- 
spect to a cold {T—0) nucleus at normal nuclear density. 
In this classical model, we adopt the similar definition 
for the ground state, ie. it corresponds to a cold nucleus 
at zero temperature and normal nuclear density where 
there is no kinetic energy and so that the ground state 

energy per nucleon is — e„^p ^'^ . Practically, N^-^ is de- 
termined by the geometry and is taken as the maximum 
bond number of unlike nucleons for disassembling source 
as it approaches to zero temperature and normal nuclear 
density as possible. 

Figure 1 shows the average multiplicitcs of the emit- 
ted neutrons, protons, charged particles and intermediate 



mass fragments, the average values of the largest frag- 
ment masses and the average excitation energies per nu- 
cleon of disassembling source as functions of tempera- 
ture and isospin. In point of temperature dependence, 
Nn, Np, Ncp and E*/A increase while Amax decreases 
with temperature as expected. However, Nimf shows a 
rise and fall with temperature as seen in previous stud- 
ies p^-|37| , which has been explained by the onset of the 
multifragmentation. In point of isospin dependence, Nn 
shows a positive correlation with / which is obviously 
reasonable due to the more neutrons for isotopes with 
larger /. Similarly, Ni„if and Amax show minor positive 
correlation with /. Reversely, Np, Ncp and E* /A seem 
to have anti-correlation with / at the same temperature 
even though the same number of protons begin with in 
these dissociation isotopic sources. Experimentally the 
similar pictures for free neutrons and protons were found 
in ^^S -I- ^'*^'^^'*Sm reactions |3^. How to understand 
these phenomena? 

Two interpretations seem to be possible. On one hand, 
symmetrical potential term in Eq.(l) will play an impor- 
tant role in controlling the emission of nucleons and clus- 
ters. With the increasing of the isospin of the sources, 
protons will feel stronger attractive potential due to the 
neighboring neutrons sharing alike increases which re- 
sults in more bound protons for disassembling sources 
with higher isospin. Viceversa, neutrons will, on one 
side, feel the stronger repulsive interaction due to more 
neutron-neutron pairs and, on other side, feel smaller at- 
tractive potential because of the decreasing of the av- 
erage assortative number of nearest neighboring protons 
for a certain neutron with the increasing of the isospin of 
sources, both reasons will lead to produce more unbound 
neutrons for disassembling sources with higher isospin. 
This explication is similar to the explanation to nucleon 
emission in isospin dependent transport model [Q . How- 
ever, another interpretation based on the excitation en- 
ergy seems to be possible, too. In a previous work we 
explored that the excitation energy can be viewed as a 
scaling quantity to control nuclear disassembly If 
we made the mapping from T to E* /A (see Fig. If) and 
re- plot the Nn, Np, Ncp, Nimf and A^ax as a function 
of excitation energy per nucleon instead of temperature 
in Figure 2. Clearly, the rule for average multiplicity of 
neutrons is not changed, ie. the higher the excitation en- 
ergy and/or isospin, the larger the Nn. But for protons, 
charged particles, intermediate mass fragments and the 
largest fragment masses, their average values for four iso- 
topes nearly merge approximately into single curves. In 
this case, the role of symmtrical potenial is transferred 
into the excitation energy of different sources and the 
isospin dependence of Np, Ncp, Nimf and A^ax at the 
same temperature showed in Fig.lb-e is nothing but the 
excitation energy dependence. 
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2. Slopes of Average Values to Temperature 

In Fig.l the average values of most physical quanti- 
ties increase or decrease monotonically with temperature 
except for the IMF. Is it possible to obtain much infor- 
mation from this figure? Figure 3 gives slopes of 
Np, Ncp, Nirnf, Araax and E* / A to temperature. Obvi- 
ously, each slope has a peak or valley around T ~ 5MeV. 
In such turning temperature, some features appear: (1) 
the emission of light particles and complex fragments in- 
creases rapidly within a narrow temperature range, re- 
flecting that the phase space is opening in the largest 
extent in that time; (2) the decrease of the largest frag- 
ment size reaches to the valley value for such a finite 
system. Physically the largest fragment is simply related 
to the order parameter pi - pg (the difference of density 
in the "liquid" and "gas" phases). In the infinite mat- 
ter, the infinite cluster exists only on the "liquid" side 
of critical point. In a finite matter, the largest cluster is 
present on both sides of critical point. In this calcula- 
tion, a valley for the slope of A^ax to temperature may 
correspond to a sudden disappearance of infinite cluster 
{"bulk liquid") near critical temperature; (3) the spe- 
cific heat C^/A (ie. d{E* / A) / dT) has a peak value for 
such a finite system. All these features are consistent to 
the concept of phase transition and critical phenomenon 
according to the thermodynamical theory. Hence, these 
slopes can be viewed as a characteristic judgement to 
critical phenomenon as other critical observables can do 
(see Sec. D). Moreover, the fact that the turning tem- 
peratures deduced from these slopes for the four isotopic 
disassembling sources nearly locate at the same temper- 
ature of ~ 5 MeV, independent of the isospin, illustrates 
that the critical temperature is insensitive to the isospin 
of the disassembling sources. This conclusion is consis- 
tent to the results deduced from other critical observables 
as shown below (see Sec. D). 

3. The Largest Fluctuations 

Furthermore, the largest fluctuations are also found 
around critical temperature in the same calculation. 
Fig. 4 illustrates that RMS width [a) for the multiplic- 
ity distributions of neutrons, protons, charged particles 
and intermediate mass fragments, for the distributions of 
the largest fragment masses and excitation energies per 
nucleon. The variances of multiplicity distributions of 
neutrons and protons have maximum values or tends to 
saturation around 5 ~ 5.5 MeV, and the variances for the 
distributions of CP, IMF, A^ax and E*/A show peaks 
at the same critical temperature. The peaks in the criti- 
cal point correspond to the cusp of the variances, ie. their 
first order derivations are discontinual across the critical 
temperature which is also an indication of phase transi- 
tion Noting that the fluctuation of A^ax is related 
to the compressibility of the system. These features are 



consistent with the critical point behavior where there 
should be the largest fluctuation in terms of statistical 
theory Similarly to the results of slopes, the deduced 
critical temperature from the largest fluctuation is al- 
most the same regardless of the isospin of the dissociation 
sources. 



B. Ratios between Two Light Particles 

Besides the multiplicity of emitted particles, the ratio 
between two light particles is probably suitable to inves- 
tigate the isospin effect. Fig. 5 shows the temperature 
dependences of the isotopic ratios: R{p/d) between the 
yield of protons to that of deuterons and R{d/t) between 
the yield of deuteron to that of triton, the isobaric ratios: 
R{t/^He) between the yield of triton and that of helium- 
3 and R{^He/^Li) between the multiplicity of helium-6 
to that of lithium-6, and the ratios R{^X/'^'^IY) be- 
tween the light particles having one proton and mass 
number difference: R{d/^He) between the multiplicity of 
deuterons and helium-3 and Rit/^He) between the mul- 
tiplicity of triton to that of helium-4 for the disassembly 
of Xe. For the isotopic ratios, R{p/d) and R{d/t), the 
particle in denominator has one more neutron than that 
in numerator, ie. the ratios could reflect the extent of 
neutron-poor of the products. Clearly, these ratios de- 
crease with the increasing of isospin as shown in Fig. 5a 
and 5b. It is consistent with the N/Z systematics of the 
disassembling source. R(p/d) and R(d/t) exists a wide 
valley around 5 MeV and then increases with temper- 
ature above 5 MeV (near the critical temperature), a 
similar result was observed experimentally in energetic 
Au -I- Au collisions These results are qualitatively 

consistent with a rise and fall in the relative production 
of IMF as shown above (see Fig. Id) . In a simple 

coalescence picture this is consistent with a lower baryon 
density at freeze-out for the most violent collision and 
disassembly [Q. For the isobaric ratios, R{t/^He) and 
R{^He/^Li), the particle in denominator has one more 
proton than that in numerator, ie. the ratio will indicate 
the extent of neutron-rich of the products. Fig. 5c and 5d 
illustrates that the larger the isospin, the higher the ra- 
tios. It has the same trend to N/Z of emitter, too. Their 
values of ratios become saturation at higher temperature. 
For the ratios of R{d/^He) and R{t/^He), the particle 
in denominator has the same neutron as that in numer- 
ator but has one more proton, ie., this ratio can show 
the extent of proton-deflcient or neutron-rich of products. 
Fig.5e and 5f show that these ratios can also imply the 
extent of isospin of the emitted source. The dependence 
of temperature is similar to R{t/^He) and Ri^He/^Li). 
Besides the above ratios, another understandable quan- 
tity is the ratio of yield of emitted neutron to protons, 
R{n/p) (not plotted here because both multiplicities of 
neutrons and protons have been shown in Fig.l), which 
increases with the isospin of the source and tends to sat- 
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uration at higher temperature. 

In order to display the isospin systematics of the above 
ratios clearly, we plot the isotopic ratios of R{p/d), 
R{d/t), R{^He/a) and R{^Li/'^ Li) as a function of the 
isospin / in a fixed temperature, eg. around critical tem- 
perature ~ 5 MeV in Fig. 6a. Obviously all the ratios 
decrease with {N — Z)/A which is consistent with the 
isospin systematics of disassembling source. 



C. Double Ratios between Two Pair Light Isotopes 

After investigating the multiplicities of LP and the ra- 
tios between one pair LP, it is interesting to study the 
double ratios between two pair LP. In the assumption of 
thermal and chemical equilibrium, S. Albergo et al. ^] 
derived a formula of temperature by the double ratios 
of isotope. The double ratios of isotope can cancel out 
chemical potential effects and offer a particular promis- 
ing technique of temperature determination p3|-|50t . But 
due to side-feeding of light particles of excited primary 
fragments and other complicate effects in experimental 
measurements, the determination of temperature from 
double ratio must be taken carefully l4^ , |5l| -|53|] . How- 
ever, here we restrict to discuss the double ratios between 
different light isotopes rather than the isotopic tempera- 
ture. 

Fig. 6b shows the double ratios between two pair iso- 
topic ratios, namely R{pd—He) — ^^^"yj^^ — , R{pd~Li) 

- ivStii' - = nZ:^M'J "^('^' - - 

N^Sl,^ ^ and R{He-Li) = ^^f^ff^- Different from 
the isotopic ratios of Fig. 6a, the isospin effect of double 
ratios is washed out as shown in Fig. 6b. The same pic- 
ture was found experimentally in central collision of ^^^Sn 
+ "^Sn and i24Sn + 1243^^ 4Q MeV/u and intepreted 
with the Expanding Emitting Source model |Q . This in- 
dependence on isospin is consistent with the attainment 
of full chemical equilibrium for each disassembling source 
at the common temperature, which probably reflects that 
Albergo's temperature is insensitive to the isospin of the 
system. 



D. Fragment Distribution and Related Critical 
Observables 

1. Fragment Distribution 

In light of previous studies on the fragment distribu- 
tion, we will use observables based on fragment to signal 
the onset of phase transition, eg. the effective power law 
parameter r, the second moment S2 of fragment distri- 
bution and the informtaion entropy. Before presenting 
these critical observables, the mass distribution of frag- 
ments are shown at T = 4, 5, 6 and 7 MeV for ^^^Xe in 
Fig. 7. Clearly the disassembly mechanism evolves with 



the nuclear temperature. A few light clusters are emit- 
ted and the big residue reserves at T = 4 MeV which 
indicates typical evaporation mechanism. With the in- 
creasing temperature, the shoulder of mass distribution 
of clusters occurs due to the competition between the 
fragmentation and the evaporation. This shoulder disap- 
pears and the mass distribution becomes power law shape 
at T = 6 MeV, corresponding to the multifragmentation 
region. When the temperature becomes much higher, the 
mass distribution becomes much steeper indicating that 
the disassembling process becomes more violent. Power 
law fits, Y(A) c>c A~'^ , for these mass distributions have 
been introduced with the lines as shown in Fig. 7. 



2. Effective Power Law Parameter 

It has already been observed that a minimum of power 
law parameter Tmin exists if the critical behavior takes 
place 56 . Fig. 8a displays r parameter as a function 
of temperature for Xe nuclei with the different isospin. 
The minimums of r parameters in Fig. 8a locate closely 
at 5.5 MeV for all the systems, which illustrates its mi- 
nor dependence on the isospin. In other words, there 
could be a universal scaling for mass distribution regard- 
less of the size of disassembling source when the critical 
phenomenon takes place. However, t parameters show 
different values outside the critical region for nuclei with 
different isospin, eg., r decreases with isospin when T > 
5.5 MeV (multifragmentation region). 



3. Campi's Second Moment 

In Fig. 8b we give the temperature dependences of 
Campi's second moment of the fragment mass distribu- 
tion 15711, which is defined as 



So - 



J2i=iA7nax Afni{Ai) 



A 



(3) 



where ni{Ai) is the number of clusters with At nucleons 
and the sum excludes the largest cluster Amaxi A is the 
mass of the system. At the percolation point S2 diverges 
in an infinite system and is at maximum in a finite sys- 
tem. Fig. 8b gives the maximums of S2 around 5.5 MeV 
for different isotopes, respectively. Again, the critical be- 
havior occurs at the same temperature, independent of 
the isospin, as r reveals. 



4- Information Entropy 

Fig.8c plots the information entropy H as a function 
of temperature for Xe isotopes. The information entropy 
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was, for the first time, introduced by H. Shannon in in- 
formation theory |5^ ]. In high energy colhsions, multi- 
particle production proceeds on the maximum stochas- 
ticity, ie. they should obey the maximum entropy prin- 
cipal. This kind of stochasticity can be quantified via 
the information entropy. In the different physical con- 
dition restrain, the information entropy can be defined 
with different stochastic variables. For examples, the in- 
formation entropy for particle multiplicity can be defined 
as 

H = -Y,pMp^)- (4) 

i 

Here pi can be defined as an event probability having 
produced particles, the sum is taken over all multi- 
plicities of products from the disassembling system. H 
reflects the capacity of the information or the extent of 
disorder. Here we introduce this entropy, for the first 
time, into the search for liquid gas phase transition and 
critical phenomenon in the nuclear disassembly. Fig. 8c 
shows that the entropy H has peaks close to 5.5 MeV 
for all isotopes. These peaks indicate that the opening 
of the phase space in the critical point is the largest. In 
other words, the system at the critical temperature has 
the largest fluctuation/stochasticity which leads to the 
largest disorder and particle production rate. Beyond 
the critical point, the information entropy H increases 
with the isospin. 

5. Caloric Curve and Specific Heat 

Another direct quantities to illustrate the phase tran- 
sition are caloric curve and specific heat. For an ideal 
infinite nuclear matter the first order liquid gas phase 
transition will reveal a temperature plateau in a wide 
range of excitation energy, which corresponds to an infi- 
nite specific heat in this constant temperature. But for a 
finite nuclear matter the plateau of caloric curve and the 
sharpness of specific heat will be largely smoothed. But 
assumably, the specific heat might keep the imprint of 
the phase transition well and probably reveal a salience 
with the temperature when phase transition and critical 
behavior holds. If exchanging the X and Y axis of Fig. If, 
the so-called caloric curves for Xe isotopes are obtained. 
For such finite nuclear systems, caloric curves present 
that temperature increases generally with excitation en- 
ergy but a temperature plateau domain is always absent 
in this lattice gas model. On the other hand, at a fixed 
temperature, the excitation energy per nucleon decreases 
with isospin of dissociation isotopic sources, which is re- 
lated to the bond numbers of unlike nucleons. Quali- 
tatively, the higher the isospin of disassembling isotopic 
source, the larger the iVj^ and then the smaller the exci- 
tation energy per nucleon in terms of Eq. (2). Due to the 
caloric curves are not linear, the slopes of excitation en- 
ergy to temperature for these curves, ie. the specific heat 



per nucleon at constant freeze-out volume (or density) 
Cy/A for Xe isotopes is useful to locate the critical point. 
Fig.Sf shows the definite peaks of C^/A exist around 5 
MeV for each disassembling sources, corresponding to the 
onset of the critical point behavior. Due to the infiucnce 
of isospin on the caloric curves, C^j/A shows a similar 
anti-correlation with isospin especially below the critical 
temperature. 

IV. SUMMARY 

In conclusion, the isospin effects of particle emission 
and critical point behavior are explored for Xe isotopes at 
a fixed moderate freeze-out density in the frame of isospin 
dependent lattice gas model. Four sets of probe are tes- 
tified to investigate such effect. First set of probe is the 
multiplicities of of LP, like Nn and iVp, and the multiplic- 
ities of charged particles N^p and of intermediate mass 
fragments Nimf, and the average mass for the largest 
fragments A^ax and the mean excitation energy per nu- 
cleon E* /A; The second one is the ratios between two LP, 
namely the isotopic ratios, like R{p/d) and R{^He/a), 
or isobaric ratios, like R{t/^He) and R{^He/^Li), or 
ratios of two light particls between which there exists 
one proton and mass number difference, like R{d/^He) 
and R{t/a); The third one is the double ratios between 
two isotope-pairs, hke R{pd — He),R{pd — Li),R{dt — 
He), R{dt — Li), and R{He — Li); The fourth one is the 
critical observables, like the effective power law param- 
eter T, the information entropy H, the Campi's second 
moment 5*2, the caloric curves and the specific heat per 
nucleon Cy/A. The calculation illustrates that Nn, Np 
and ratios of light particles show the strong dependences 
on the isospin of the dissociation source, but it's not the 
case for the double ratios of light isotopes. The former 
reflects directly the chemical composition of the source 
and the latter indicates the chemical equilibrium of the 
source, which is not self-contradict but reflects features 
of the sources on different sides. Meanwhile, the criti- 
cal temperature for a chain isotopes determined by the 
extreme values of t, H and 5*2 is also insensitive to the 
isospin of sources. This conclusion is not contradict with 
the previous studies on the isospin dependence of critical 
temperature like in |2^] , where the span of isospin is from 
symmetrical nuclear matter to pure neutron matter. If 
we only look a small span of isospin for the experimen- 
tally measurable medium size isotopes, like ^^^~^^^Xe, 
the change of critical temperature is neglectable. This 
conclusion might indicate that it will be difficult to search 
for the isospin dependence of critical temperature which 
signals the liquid gas phase transition for medium size 
nuclei in the experimental point of view. In addition, 
values of power law parameter of cluster mass distribu- 
tion, mean multiplicity of intermediate mass fragments 
(IMF), information entropy (H) and Campi's second 
moment (S'2) also show minor dependence on the isospin 
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of Xc isotopes at the critical point. In contrary, some 
isospin dependences of the values of t, _ff and 5*2 will 
reveal outside the critical region. Noting that the infor- 
mation entropy is, for the first time, introduced into the 
search for the liquid gas phase transition in this work. Be- 
sides, the slopes of A^„, Np, Ncp, Nimf, A^ax and E*/A 
to temperature and the variances of the distributions of 
these quantities are explored as additional judgements 
for critical phenomenon. Similaring to the insensitivity 
of critical temperature to the isospin from the extreme 
values of r, H and ^2, the critical temperature deduced 
from the above slopes and variances also shows indepen- 
dence on the isospin of the disassembling sources. It will 
be interesting to confront these conclusions with the fu- 
ture experiments. 
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Fig.7: Mass distribution of ^^SXc at T = 4, 5, 6 and 7 
MeV. The linos are the power-law fit. 

Fig. 8: Critical obscrvables: t parameter from the 
power law fit to mass distribution (a), Campi's second 
moment (b) and information entropy (c) as functions of 
temperature and isospin of disassembling sources. 



Fig.l: Average multiplicities for emitted neutrons (a), 
protons (b), charged particles (c), intermediate mass 
fragments (d), average values of the largest fragment 
masses (e) and of excitation energies per nucleon (f) as 
functions of temperature and isospin of Xe. 

Fig. 2: The same! as Fig.l but as a function of excitation 
energy instead of temperature. 

Fig.3: Slopes of 7V„ (a), Np (b), N^p (c), iV„„/ (d) and 
■^max (e) to temperature. The symbol is the same as in 
Fig.l 

Fig. 4: Variances for the distributions of iV„ (a), Np 
(b), N^p (c), Ni^f (d), 

-^max (c) and E* / A (f) versus 

temperature. 

Fig. 5: Isotopic ratios R{p/d) (a) and R{d/t) (b), iso- 
baric ratios R(t/^He) (c) and R(^He/^Li) (d), and ratios 
R{d/^He) (e) and R{t/^He) (f) as functions of temper- 
ature and isospin of Xe. 

Fig.6: lotopic ratios of R{p/d), R{d/t), R{^He/a) and 
R{^ Li Li) (a) and double ratios R(pd—He), R(pd—Li), 
R{dt — He), R{dt — Li) and R{He — Li) (b) as a function 
of the isospin {N — Z)/A around critical temperature 
~ 5 MeV. The symbols are illustrated in (a) and (b), 
respectively. 
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